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It is shown that a strong impurity potential induces short-range antiferromagnetic (ferrimag- 
netic) order around itself in a Hubbard model on a half-filled honeycomb lattice. This implies 
that short-range magnetic order is induced in monolayer graphene by a nonmagnetic defect 
such as a vacancy with full hydrogen termination or a chemisorption defect. 
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Physical properties of graphene (a monolayer graphite 
sheet) have attracted much interest since it was fab- 
ricated, 1 magneto-transport properties were studied, 2,3 
and various peculiar properties such as the half-integer 
quantum Hall effect 4-6 were found. Although the main 
interest is in its transport properties at present, mag- 
netic properties are also worth studying. Indeed, inter- 
est in magnetism in carbon-based materials has recently 
surged because of not only its fundamental importance, 
but also possibility of application in new technologies. 7 

Electrons in graphene can be described by a tight- 
binding model on a honeycomb lattice. The purpose of 
this paper is to show that short-range antiferromagnetic 
(AF) (more precisely, fcrrimagnctic) order is induced 
around an impurity introduced on a half-filled honey- 
comb lattice. 

Magnetism in graphite systems has been studied in 
many experiments, 8 ' 9 although the issue has been con- 
troversial; the effect of magnetic impurities cannot be 
completely discarded. Apart from magnetic impurities, 
the defect-induced mechanism is the most probable 
mechanism of magnetism in carbon-based materials. For 
example, a free dangling bond associated with a vacancy 
will have a net magnetic moment, which may induce a 
magnetic state. Even if a defect itself is nonmagnetic, it 
can cause spin polarization of 7r electrons around it. It has 
been known that localized edge states are generated at a 
zigzag edge and they can lead to a magnetic state. 10-12 
In previous studies, 13 ' 14 we showed that a strong impu- 
rity potential in a half-filled honeycomb lattice induces 
a localized state (zero mode) around it 15-17 and that it 
causes enhancement of staggered susceptibility. We did 
not find any (short-range) order, however, because we did 
not consider electron-electron interaction. In this study, 
considering electron-electron interaction, we show that 
short-range magnetic order is indeed induced around a 
strong-impurity site in a half-filled honeycomb lattice. 

We start from a tight-binding model on a honeycomb 
lattice. We consider only the nearest neighbor transfer t, 
and then the hamiltonian is given by 



Ho 



(nA.n'B) 



T CnBo 



h.C.], 



(1) 



where c naa [Cnaa) i s an annihilation (creation) opera- 
tor of an electron with the spin a on sublattice a in 
the nth unit cell, and (nA, n'B) stands for the pairs of 



the nearest sites. We denote the number of the unit cell 
by N u = L 2 . The number of the lattice points is JVl 
2 Y . = 2L 2 . In addition to Ho, we consider electron- 
electron interaction. We consider the on-site Coulomb 
interaction U. Long-range Coulomb interaction is poorly 
screened in graphene. However, because it must be the 
short-range part of the Coulomb interaction that plays 
the most important role in inducing magnetism, we con- 
sider only the on-site Coulomb interaction and use the 
following interaction hamiltonian H', 



H' = U ^ n na t fl na L, 
n,a 



(2) 



where h nacT = c na a c naa- We thus consider a Hubbard 
model on the honeycomb lattice. 

Furthermore, we introduce impurities. We consider a 
short-range impurity potential. If an impurity is on sub- 
lattice a in the nth unit cell, its effect is expressed by 



V(n,a) 



naa^naa ■ 



(3) 



The limit of u — > oo represents a strong impurity such as 
a vacancy or a hydrogen chemisorption defect. 18 In that 
case, the total number of the lattice points is N-^ = 2 AT — 
Aj, where N\ is the number of impurities. We consider 
the half-filled case, i.e., the case where the number N c 
of electrons is equal to A L , in this study. Actually, we 
consider the case with Ni = 2. 

We resort to a mean field approximation. The interac- 
tion term is approximated as 



n.a 



U ^(flnal) (final) + U ^ (^na-a)fl 
n.a. a 



(4) 



The mean field hamiltonian H m f is then given by, apart 
from a c-number term, 



Tin 



naa 



-naa ' h naa 



Ho + V, 



(5) 



where e naa = V (fi n a-a) ■ By diagonalizing H m f, we self- 
consistently determine 2Al parameters, e nQff , or, cquiv- 
alcntly, (h naa ), from which we can readily calculate the 
electron density n na = (final + f-nai) and the spin den- 
sity nina = (final ~ "rial) at each lattice point. 

Before showing the results in the presence of impuri- 
ties, we briefly study the impurity-free case. We define 
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static spin susceptibilities X±(?) by 



1 f 13 

X±W = -TT / dr ([m A {q,T) ±m B (q 1 T)} 

iV u Jo 

x [m A {-q,0)±m B (-q,0)] ), (6) 



where 

X+(0) is uniform susceptibility and X-(0) is staggered 
susceptibility. In Fig. 1, we show the momentum depen- 
dence of X±(q) f° r U = at absolute zero, T = 0. It can 
be seen that staggered susceptibility \-(q = 0) takes the 
maximum value. 

In the random phase approximation (RPA), spin sus- 
ceptibilities can be written as 

xL 0) ( 9 ) 



x±(q) 



i-fxi 0) («y 



(8) 



from which we can see that the paramagnetic state is 
unstable towards an AF state at U > Uc? = 2 /%_ (0) ~ 
2.23*. 
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Fig. 1. Momentum dependence of spin susceptibilities x+(<?) 

(dotted curve) and X— (<?) (solid curve) of the half-filled hon- 
eycomb lattice at U = and T = 0. 



(ni,ni) and and rii = n\ + (L/2, L/2), and impose the 
periodic boundary condition. As can be seen from Fig. 
2, the (shortest) distance between the two impurity sites 
is L/2. 

First, we study the case with u — ► oo, i.e., JVl = 
2N U — 2, and deal with the case where both impurities 
are on the same sublattice, say sublattice A; we denote 
this case by (A, A). In this case, the ground state is four- 
fold degenerate at U = 0, because four degenerate zero 
modes 4>^{r){a =T, j and k = 1,2) appear just at the 
chemical potential ijl (see Fig. 3 (a)). At a finite U, two of 
the zero modes are occupied by equal spins in the ground 
state (Hund's rule), that is, the ground state always has a 
total moment M t = m A + m B = 2 where m a is the total 
moment of the sublattice a, m a = ^2 n (n na -\ — n> na i). 
This conclusion is in harmony with Lieb's theorem. 19 
Because the zero-mode wave function (jy^J (r) is localized 
around impurity sites, the spin density is also localized 
around impurity sites. In Fig. 4, we show the spin density 
around an impurity site. It can be seen that short-range 
fcrrimagnctic order is induced around the impurity site 
and that the spin density on sublattice B is larger than 
that on sublattice A, the sublattice where the impurity 
is. 20 As a finite total moment M t implies, the direction of 
the net spin density around one impurity site is parallel 
with that around the other impurity site. This is eas- 
ily understood as a manifestation of the dominant AF 
tendency of the impurity-free half-filled honeycomb lat- 
tice. The net majority moments around both impurities 
are mainly on sublattice B, and the nearest-neighbor AF 
interaction implies a ferromagnetic interaction between 
spins on the same sublattice. We show the dependence 
of the moment at each lattice point on the distance from 
impurity sites in Fig. 5 (a). It can indeed be seen that the 
moment is localized around impurities and ferrimagnetic 
short-range order develops. 




Fig. 2. Honeycomb lattice of size L X L used in calculations. The 
periodic boundary condition is imposed. Solid dots represents 
the positions of the unit cells where impurities are placed, and 
open dots are their images, (i = 1, 2) is the unit vector of 
length a, where a is the lattice constant. The dotted line shows a 
unit cell which includes an A sublattice point and a B sublattice 
point. 



Now, we study the impurity effect using a finite size 
lattice. The lattice with L x L unit cells is shown in Fig. 
2. We place two impurities in the unit cells at ri\ = 



(a) (A,A) 
U=0 T 



U>0 
(b) (A,B) 



Fig. 3. The ground state spin configuration in the presence of two 
vacancies (a) on one of the sublattices and (b) on sublattices A 
and B. (a) At U = 0, four zero modes (two per each spin) appear 
at the chemical potential fi. At U > 0, two of them are occupied 
by equal spins (up spins in this figure), (b) Zero modes appear 
symmetrically around the chemical potential fj,, and two of them 
are occupied by an up-spin and a down-spin. 



The dependence of the square of staggered moment 
M s = m A — m B on U is shown in Fig. 6; M s smoothly 
changes from 2 as U increases from zero. Total sublattice 
moments tua and m B are given by m A = (2 + M s )/2 and 
mB = (2 — M s )/2 (M s < 0). This means that no moment 
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is induced on sublattice A as U — > 0. This is because the 
zero-mode wave function vanishes on sublattice A. 14 ' 15 
The sublattice moment on sublattice A is induced by 
a finite U . This can be also seen in Figs. 4 (b) and 5 
(a). The local moment M^{i max ), defined by M[(£ max ) = 
^2n ( m nA — mnB) where the summation is restricted to 
those lattice points within the distance of f max a from an 
impurity site, is also shown as a function of U in Fig. 7. In 
contrast to the staggered moment M s , the local moment 
varies only modestly. 21 It is the staggered moment that 
is enhanced by the interaction. 




(a)U/Ucr = 0.5 3 



(b)U/Ucr = 0.09 





Fig. 6. The dependence of the square of staggered moment M B on 
the interaction strength U; the linear size of the lattice is L = 40 
(solid curves) and L = 64 (dotted curves). M s varies smoothly 
in the case (A, A) where both impurities are on sublattice A, 
while it evolves from zero at an L-dependent critical U cr (L) in 
the case (A,B) where one impurity is on sublattice A and the 
other is on sublattice B. The upper inset shows the dependence 
of M s on U in a wider range for L = 40. The lower inset shows 
the dependence of U C r(L) on 1/L. The critical value Uc? for the 



bulk honeycomb lattice is U t 



(0) 



2.23t. 



Fig. 4. Spin density around an impurity site at (a) U = 1.2t = 
0.53t/ c ( r 0) and (b) U = 0.2t = 0.09u£ ] for L = 50. Area of each 
dot is proportional to the magnitude of the spin density at each 
lattice point. Open dots represent positive magnetic moments 
and filled dots negative magnetic moments. 




Fig. 5. Moments — rra„'A (solid dots) and m n ig (open dots) in 
the unit cells at n! = n\ + (n, n) at U = l.Ot for L = 64. (a) Im- 
purities are on the A sublattice points at n± and at n\ + (32, 32). 
For comparison, —m n /j^ at U = 0.5t is also shown (stars), (b) 
Impurities are on the A sublattice point at n\ and on the B 
sublattice point at n\ + (32,32). 



Next, we study the case where one vacancy is on sub- 
lattice A, and the other is on sublattice B; we denote this 
case by (A,B). In this case, the zero modes appear pair- 
wise above and below the chemical potential (see Fig. 
3 (b)). At any U, two of the zero modes arc occupied 
by up- and down-spins and the total moment M t always 
vanishes in harmony with Lieb's theorem. 19 At U = 0, 
the occupied zero-mode wave function </>g~^(r) takes a fi- 
nite value only on sublattice a around an impurity site on 
sublattice a. 14 ' 22 The spin density at each site vanishes, 



1 




Fig. 7. Local spin moment M t (£ max ) as a function of U for L = 
40: £ ma x = 16a (solid curves), 8a (dotted curves), and 4a (dashed 
curves). 



because (f> ^ (r) — (p^7'(r). As U increases, this identity 
is violated. For example, (r) | takes a larger value on 

sublattice B than ( r )l> ano ^ vice versa 011 sublattice 
A. This results in a net up spin density around the im- 
purity site on sublattice A and a net down spin density 
around the impurity site on sublattice B. The spin den- 
sity around an impurity site at U = 1.2t is very similar 
to the one in the case (A, A) at the same value of U, i.e., 
Fig. 4 (a). We again find short-range fcrrimagnctic order 
around a vacancy site. See also Fig. 5 (b). In this case, the 
induced net moment around the impurity site on sublat- 
tice A is antiparallcl to that around the impurity site on 
sublattice B. This is also understood as a result of the in- 
herent AF tendency of the half-filled honeycomb lattice. 
Actually, the net moment induced around one impurity 
is equal in magnitude with that induced around the other 
impurity. 

We find that, for a finite L, the linear size of the lat- 
tice, there is a critical value of U, U CT (L), above which 
staggered moment takes a finite value, as can be seen 
from Fig. 6. Near U CT (L), the dependence of M s is well 
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approximated as M s oc yi7 — U C1 -(L). The dependence 
of U CT (L) on L is shown in one of the insets of Fig. 
6. The previous results 13,14 strongly suggested that the 
staggered susceptibility on the non-interacting half-filled 
honeycomb lattice diverges in the presence of a strong 
impurity potential. This implies that U cr (L) vanishes as 
L — ► oo, that is, a finite staggered moment is induced 
by an infinitcsimally small value of U . Moreover, as can 
be seen from Fig. 6, the \M S \ — U curve in the (A.B) 
case merges that in the (A, A) case as U increases, and 
it merges at a smaller U at a larger L. This also implies 
that U CI (L) vanishes as L — > oo. Indeed, it will make 
no difference whether two impurities are on the same 
sublattice or on the different sublattices, if the distance 
between them goes to infinity. From these arguments, 23 
it is clear that U ct {L — > oo) vanishes, although we can- 
not definitely confirm U CI (L — > oo) = from the present 
calculation because of the limited size of the used lattice. 
The dependence of the local moment M*(£ max ) on U is 
shown in Fig. 7. It also merges Mjf (£ max ) in the (A, A) 
case as U increases. 

The local moment formation is closely related to the 
presence of the zero modes, i.e., it is specific to the half- 
filled honeycomb lattice. By doing similar calculations 
for an off-half-fillcd honeycomb lattice and a half-filled 
square lattice (without nesting), we indeed find that no 
local moment formation occurs at U much smaller than 
Uci . Vacancies on the same sublattice points in a row 
generate a zigzag edge (and a bearded edge) where fer- 
rimagnetc order is realized. 10, 12 This order at a zigzag 
edge is thus closely connected to the short-range order 
around a vacancy studied here. 
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Fig. 8. Region where local moments are formed on the u — U 
plane in the (A,B) case obtained for L = 20 (open dots) and 50 
(filled dots). Electron number N e is equal to — AT; = — 2. 
At U > u£' , the AF state is realized. 



Finally, we briefly discuss the case with finite u. In Fig. 
8, we show the region where a local moment is formed 
around an impurity site in the case (A,B). At «<T0i, 
results converge to those at u — ► oo. At smaller u, the 
region where a local moment is formed is rather limited. 
Finite size effects are clearly seen, and the detailed study 
of them is left for a future study. 

Now, we discuss the applicability of the present result 
to actual graphene. We do not know the precise value 
of U appropriate for graphene, and it is also difficult to 
estimate the value of u that represents the effect of a de- 
fect. Moreover, we have neglected the effect of long-range 
part of the Coulomb interaction and the possible local 



distortion around a defect. Yazyev and Helm recently 
found fcrrimagnetic short-range order around a vacancy 
and also around a hydrogen chemisorption defect using 
a first-principles calculation. 24 This implies that actual 
graphene with nonmagnetic defects is in the region of 
Fig. 8 where local moments develop around impurities 
and that the short-range order is caused by the same 
mechanism discussed in this work. 

We find that the energy gain per vacancy due to the 
moment formation is approximately O.Olt, 25 i.e., 300 K 
in graphene. This implies that it is difficult to observe 
local moment formation around a single vacancy at room 
temperatures, and that it can be observed only at low 
temperatures. It is possible that a moment develops more 
stably near a few defects close to each other. This is also 
a subject of a future study. 

A local magnetic moment can also be formed around 
a defect in other carbon-based materials. For example, 
oscillation of spin density was found around a vacancy in 
a (10,0) carbon nanotube with a first-principles calcula- 
tion. 26 This may also be caused by the same mechanism 
discussed here. 

To conclude, we have studied the effect of nonmagnetic 
defects on magnetism of electrons on the half-filled hon- 
eycomb lattice and found that short-range ferrimagnctic 
order is induced and a local moment is formed around 
a defect. This strongly implies that a nonmagnetic de- 
fect in graphene indeed induces a local magnetic moment 
around it. 
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